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The low-energy fluctuations of a trapped, interacting quasi one-dimensional Bose gas are studied. 
Our considerations apply to experiments with highly anisotropic traps. We show that under suitable 
experimental conditions the system can be described as a Luttinger liquid. This implies that the 
correlation function of the bosons decays algebraically preventing Bose-Einstein condensation. At 
significantly lower temperatures a finite size gap destroys the Luttinger liquid picture and Bose- 
Einstein condensation is again possible. 

PACS: 03.75.Fi, 05.30.jp, 32.80.Pj, 67.90.-fz, 71.10.Pm 
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The experimental realization of Bose-Einstein conden- 
sation (BEC) in atomic vapors of ^'''Rb and ^"^Na 
1^,^ has attracted a lot of interest Q. Recently, a 
highly anisotropic, quasi one-dimensional trap has been 
designed [||. Up to now, the possibility of BEC in 
one dimension has mainly been discussed for the non- 
interacting Bose gas The role of dimensionality 
has been carefully examined for the ideal bose gas by van 
Druten and Ketterle j|] . In one dimension the interaction 
between bosons plays an essential role due to the strong 
constraint in phase space The question of BEC in 
a quasi one-dimensional system is therefore more com- 
plicated. The purpose of this paper is to demonstrate 
that under suitable experimental conditions the low en- 
ergy excitations of this system are described by a Lut- 
tinger liquid (LL) |0 model. The superfluid correlations 
of a LL decay algebraically and the system is not Bose 
condensed. At much lower temperatures which are de- 
termined by the extension of the trap in the longitudinal 
direction the spectrum of the phase fluctuations is again 
cut off by finite size effects and the bosons could condense 
again. 

The realization of a Luttinger liquid in a one- 
dimensional Bose gas would be a highly non-trivial exam- 
ple of an interacting quantum liquid. Fermionic systems 
which are believed to be described by a Luttinger liquid 
include quasi one-dimensional organic metals pl] |, mag- 
netic chain compounds, quantum wires and edge states in 
the Quantum Hall Effect. While these systems are always 
embedded in a three-dimensional matrix and thus show 
a crossover to a three-dimensional behavior at low tem- 
peratures, the trapped one-dimensional Bose gas would 
provide a clean testing ground for the concept of a Lut- 
tinger liquid. 

The paper is organized as follows: First we discuss the 



circumstances under which a trapped Bose gas can be 
considered as a one-dimensional quantum system. Next 
we demonstrate in an explicit calculation that there is a 
gapless mode with a linear dispersion. We show that the 
Hamiltonian of the low-lying excitations can be identified 
as that of a Luttinger liquid and therefore the density- 
density correlation function decays algebraically. In the 
reminder of the paper we discuss the implication of the 
algebraic decay of the particle-particle correlation func- 
tion for BEC and review the properties of a Luttinger 
liquid. 

We consider the Bose gas in a cylindrical symmetric 
trap confined to the z-axis by a tight trapping potential 
in the xy-plane. If the extension L of the trap in z- 
direction is much larger than its radius R, it is justified 
to approximate the potential in the longitudinal direction 
by zero. One-dimensional physics will be dominant, if the 
temperature is much lower than the energy of the lowest 
radial excitation. The energy scale is set by huj±, with 
Lu±_ being the trap frequency 
for one-dimensionality is 
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Thus the condition 
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where T the temperature of the Bose gas. A typical value 
for uj± which has been realized in the experiments per- 
formed at the MIT by the Ketterle group (| is 27r-240i7z. 
In order to realize a one-dimensional Bose gas for this 
value of uj±_ the temperature has to be lower than l.SnK . 
Another possibility is to increase the value luj_ which 
might be more feasible experimentally. For instance per- 
manent magnets can be used to increase trap frequencies 
by more than an order of magnitude [l^ . 

Assuming that this condition for uj± can be met ex- 
perimentally, we can model the system by the following 
Hamiltonian: 
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A + C/(f)-/iU(f) (2) 



+ -J(fr (fir' i}\r) ^\r') g6{r- r") ^{f") ^P{f) , 

where m is the atomic mass, ^ is the chemical poten- 
tial fixed by the particle number N = J d^r\ip{r)\'^ and 
g = Anh'^a/m is the coupling constant, with a being the 
s-wave scattering length. We only consider repulsive in- 
teractions. U = ^muj'^{x'^ +y'^) is the trapping potential. 
The field operators ■0t(r) and tp{r) are bosonic creation 
and destruction operators. 

We now illustrate that a gapless mode for a Hamilto- 
nian like (||) exists The dynamics of tlj{r,t) is 

governed by the equation of motion 

^2 

ihdtip^~ — At{j + {U-n)ip + gij''ipij, (3) 
2m 

a one-dimensional non-linear Schrodinger equation, Eq. 

which well understood We merely illustrate in 
the following its application to the problem of trapped 
bosons. For a macroscopically occupied ground state, 
the operator ip can be considered as a classical complex 
field. Then Eq. (^) becomes the Gross-Pitaevskii equa- 
tion. We describe the complex field ip{r, t) by its density- 
phase representation: ^/j{f,t) — yj p{r, t) exp{i6{r, t)). A 
saddle point solution to Eq. (j^) is given by a constant 
phase and static density p{f,t) = Po(^) which only de- 
pends on the radius r, due to the axial symmetry of the 
problem. The solution of the Gross-Pitaevskii equation 
in a cylindrical trap and its fluctuations in the Thomas- 
Fermi approximation has been discussed in detail by E. 
Zaremba We only repeat the steps of the calcula- 

tion necessary for our arguments. Expanding in small 
fluctuations of the phase, S0, and density, 6p, around the 
saddle point solution: 

we obtain the linearized equations of motion for 6p and 
66 

ndt6e^gSp-^-\/(po\/^) , (4) 
4m po V PoJ 
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hdtSp^ —vipovse) 
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These equations possess a trivial solution {6p ^ , 
SO — const.) whose energy vanishes. This is the Gold- 
stone mode corresponding to global rotations of the con- 
densate's phase. Radial fluctuations can be ignored, be- 
cause their energy scale is set by huj±, the trap frequency. 
Thus it is justified to consider the one-dimensional limit 
where the equations simplify to 



t-2 1 

hdtSe^gSp^- d^J 

4m Po 




FIG. 1. Qualitative behavior of the low energy excitation 
spectrum (lower curve) for a one-dimensional trapped Bose 
gas. The first radial excitation (upper curve) with energy 
uj ~ uj± is also shown. The discussion in the text focusses on 
the role of the lower branch. 

The solutions are plane waves {Sp , 66 oc e'*^'^~'^'^) with 
frequencies 
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where the sound velocity is given by Vg = \/ gpa/fn which 
is the Bogoliubov value for a homogeneous Bose gas. It 
has been observed experimentally in anisotropic 3D traps 

We draw two important conclusions from this relation. 
The g^-term cannot be treated as a small perturbation 
on the energy of a non-interacting Bose gas, hence the 
smallest interaction changes the excitation spectrum fun- 
damentally. The existence of a collective mode with lin- 
ear dispersion for small g is a direct consequence of the 
interaction between particles. Only for a vanishing cou- 
pling constant g, the spectrum reduces to that of free 
particles, regardless of the ground state occupation. In 
a one dimensional trap (L 3> R) the phase-fiuctuations 
of the boson wave function destroy superfluid order due 
to phase space constraints [Q. The finite size gap in 
three dimensional traps(L « R) introduces a cut-off in 
the phase space integrals, the phase space argument does 
not apply ||l2|,|l9| and BEG is possible. As all trapped 
Bose gases are of finite size, in principle the phonon spec- 
trum remains discrete. The level splitting is only relevant 
in the limit 
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fc^r ^ hv. 
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hdt5p = — po9f 
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For the MIT trap g with the length L = 0.5 mm, this 
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temperature turns out to be roughly 10~^^K (assuming 
Na atoms) . Only in this limit the system can be in a Bose 
condensed phase. If the length L is not macroscopic the 
gap in the lowest mode will be appreciable and there is 
Bose-Einstein condensation for a finite number of parti- 
cles as pointed out by Ho and Ma iQ. We stress that 
this is due to the smallness of L and not a generic fea- 
ture of the system. A setup with L comparable to R 
is really three-dimensional. One can check that the gap 
energy for a one-dimensional trap found by Ho and Ma 
scales with the inverse axial extension of the system and 
hence disappears for large systems. We conclude that for 
L satisfying condition (|9|) there is a gapless sound mode 
which inhibits the formation of a condensate at all fi- 
nite temperatures. Nonetheless the decay of coherence is 
only weak. This is due to the fact that the system can 
be described as a Luttinger liquid as will be shown now. 

With the same approximation as for the equations of 
motion, the Hamiltonian in the long-wavelength limit is: 



H 



dz 



(10) 



where p is the number of particles per unit length and k 
is the compressibility. 

The Hamiltonian, Eq. (px|), is known as the Luttinger 
liquid Hamiltonian (l^,^^ . This concept has been mostly 
used to investigate the properties of fermionic systems in 
one dimension. The Luttinger liquid Hamiltonian, Eq. 
(p^), can be diagonalized by a Bogoliubov transforma- 
tion in terms of new bosonic creation and destruction op- 
erators 6j , bq for the long- wavelength density- fluctuation 
modes. This is possible due to the linear dispersion rela- 
tion. The Bogoliubov-transformation is given by: 



(11) 
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The b^jbq satisfy the usual boson commutation relation 

[6^,6^/] — 6q^q' and 66 and 6p form a pair of conjugate 
operators: 



[6d{z),Sp{z')]^iS{z~ z') . 
This condition fixes the functions fq and gq: 



fq - Vk|e"' 
9q = sgn{q)^ 



(13) 
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where is the parameter of the Bogoliubov- 

Transformation. Inserting the representations ( pj] ) and 
( p^ ) for 5p and 69 respectively in the Hamiltonian for the 
fluctuations the LL in terms of the new bosonic operators 
is given by: 
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with the choice exp(2ag) — h\j~pjmg. The phonon fre- 
quency is given by tjg = Vs where Vs is the sound 
velocity. 

One of the striking features of a Luttinger liquid is 
that the model has only two microscopic parameters, 
the sound velocity Vs and the compressibility k. An- 
other important property is that the correlation functions 
of the original boson operators decay algebraically in a 
Luttinger liquid. The asymptotic behavior for large dis- 
tances, z — > cx), of the boson-boson and density-density 
correlation function is given by [|lO[ : 



(*^(z)«'(0)) - 1/z 
(p(z)p(O)) - ^ r^lz 



l/r, 



(17) 
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where ry is the correlation exponent. A useful naive es- 
timate for ?/, assuming that the compressibility k ^ gp^ 
is: 



(19) 



where I b = \/hJuj\m is the magnetic length of the trap 
perpendicular to the z-axis and a is the scattering length 
of the trapped atoms. Because the interaction is weak 
we do not expect the exponent 77 to be renormalized sub- 
stantially. For current traps the exponent 77 is of the 
order ry ^ 1000 demonstrating that the phase coherence 
of the bosons decays only very weakly and is experimen- 
tally undistinguishable from true BEG |2^]. However for 
steeper magnetic traps, wj, ~ 50 kHz, particle densi- 
ties of p '--^ 10'* particles/cm and assuming a scattering 
length of 110 as for Rb the exponent 77 is 77 ~ 4 
and it should be possible to observe LL behavior. Be- 
low T ^ 0.4 nK only the linear mode is excited and the 
physics is described by LL physics. At still lower temper- 
atures, T ~ 10~^^ K the finite size gap comes into play 
]19|- 

Next we compare our results to the "two-step conden- 
sation" picture put forward by van Druten and Ketterle 
The authors consider an ideal Bose gas in a highly 
anisotropic trap. In the non-interacting system there is 
no fundamental difference between the one and three di- 
mensions except in the density of states. As soon as 
interactions have to be considered the situation changes 
drastically. Basically we have developed a more precise 
physical picture of the regime which van Druten and Ket- 
terle call the "two-step BEG" [||. Our claim is that in 
this regime the ground state is described by a Luttinger 
liquid and not by an ideal Bose gas. 

Since the Luttinger liquid model has a harmonic 
Hamiltonian, Eq. (p^), for the phase and density fluctu- 
ations, any expectation value and dynamical correlation 
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function of the boson operators in the long-wavelength 
limit can be evaluated. Luttinger liquids are well under- 
stood and many results can be carried over to the one- 
dimensional trapped Bose gas. At larger densities the 
parameters of the Luttinger liquid model will be renor- 
malized from the saddle point values by short range fluc- 
tuations and also by the three-dimensional density profile 
of the trapped Bose gas. The renormalized parameters 
can be obtained by considering more realistic interactions 
in one dimension. For a repulsive delta-function poten- 
tial the sound velocity and the compressibility have been 
obtained exactly . We are currently working on mod- 
els with longer range interactions which can be treated 
by the Density-Matrix-Renormalization-Group method. 
Results of this work will be presented elsewhere [|4| . In 
a complementary approach, we calculate the finite size 
effects in the experimental setup on the dynamics of the 
bosons 1^^. Another interesting problem which we are 
currently investigating is the response of the system to 
an impurity atom. The finite mass leads to an unusual 
behavior of the mobility ||2^,0 . Also the transport prop- 
erties should differ significantly from the conventional 
Bose condensate if the Bose gas is in the Luttinger liquid 
regime. 

To summarize, we have shown under which experimen- 
tal conditions a trapped quasi-one dimensional system 
of interacting Bosons is described by a Luttinger liq- 
uid Hamiltonian. An experimental realization of such 
a system would provide a clean laboratory for testing 
the properties of a Luttinger liquid. Its behavior devi- 
ates significantly from the noninteracting Bose gas. Un- 
like other systems which are realizations of a Luttinger 
liquid, three-dimensional effects become less important 
for lower temperatures. Moreover, it would be possible 
to tune important parameters like the density and the 
length, i.e. the trap frequency lo±^ which is impossible in 
a solid. 
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